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Abstract. The singular nature of contact problems, such as walking, makes 
them difficult to analyze mathematically. In this paper we will "regularize" the 
contact problem of walking by approximating the ground with a smooth repul- 
sive potential energy and a smooth dissipative friction force. Using this model 
we are able to prove the existence of a limit cycle for a periodically perturbed 
system which consists of three masses connected by springs. In particular, this 
limit cycle exists in a symmetry reduced phase. In the unreduced phase space, 
the motion of the masses resembles walking. 



1. Introduction 

Consider the physics of walking. The number of degrees of freedom in the system 
seem to change, depending on whether one foot is on the ground or in the air. The 
on/off nature of contact makes the analysis of such a system difficult. In this 
paper, we will "regularize" the contact problem by using smooth viscous friction 
and potential energies to model the ground, resulting in a dissipative system on a 
phase space TQ. We will then identify an M symmetry of the system and find a 
reduced order model on the quotient space TQ/M. Upon finding a hyperbolically 
stable fixed point in TQ/M., we will create a limit cycle by adding a periodic 
perturbation. This limit cycle in TQ/M. corresponds to a "relative limit cycle," or 
walking, in the original phase space TQ (see figure [T]). 
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Figure 1. This commutative diagram characterizes walking as a 
relative limit cycle in a dissipative Lagrangian system. 
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1.1. Main Contributions. The main contributions of this paper are: 

(1) We present a model of a walker in IR 2 where the no-slip condition and the 
no penetration condition have been "smeared" over a region around y — 0. 
This model will be a dissipative Lagrangian system on a tangent bundle 
TQ where Q = IR 6 . 

(2) We will determine an M-symmetry for the system which corresponds to 
translation along the x-axis. We will also state the reduced equations of 
motion on TQ/M. 

(3) We prove the existence of a stable point for our 2d walker in TQ/M.. 

(4) We will periodically perturb the system by oscillating the equilibrium 
length of the legs. This will allow us to assert the existence of a limit 
cycle in TQ/M.. 

(5) The limit cycle in TQ/M will be related to trajectories in TQ via a phase 
reconstruction formula. We will find that when the phase is non-zero, the 
trajectories in TQ correspond to "walking" as defined in this paper. 

1.2. Background Sc Motivation. This "regularized" approach is in contrast to 
the "hybrid systems" approach, where transitions between different types of phase 
spaces are given by various transition maps. The hybrid systems approach is ar- 
guably more accurate, and has yielded a number of insights and useful models. 
For example, a hybrid systems formulation was introduced in [12] where the tran- 
sition maps led to a dimension reduction; it was suggested that a limit cycle was 
approached passively. Since the work of [12] . the notion of walking as a limit cycle 
has become more common, and more sophisticated analysis have lent further sup- 
port to this idea [5] [6] . The most compelling arguments are the original videos of 
McGreer which accompany [T2] . 

On the biological side, the notion of "central pattern generators" (CPGs) has 
been hypothesized to be a fundamental neural mechanism used in walking [7J. 
These CPGs are non-localized collections of neurons which produce rhythmic ac- 
tivity, and respond to various inputs which modulate these rhythms. Therefore the 
link between CPGs and limit cycle walking is one which links periodic activation 
of the controls to periodic motion of the body. This link is used in the creation of 
simple models which can be feasibly analyzed (see for example [6]). 

Lastly, viewing walking as a limit cycle allows for great reductions in complexity. 
In particular, under weak assumptions, the existence of limit cycles in hybrid 
systems implies the existence of a reduced order model for the system as a whole 
[2J. The most recent paper, [2J, dealt with the singular nature of hybrid systems 
by using relatively weak assumptions on the transition maps to "smooth" the 
dynamics across the transition regions and prove the existence of a lower-order 
hybrid model. Our paper can be seen as a "dual approach" to [2J in that we 
regularize the transitions maps themselves. To do this we will use viscous friction 
force to approximate the no-slip condition as described in [11] and §1.6 of pQ. 
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Using our smooth system we will model walking as a limit cycle. It seems 
fair to ask "in what space does this limit cycle exist". In this paper, we will 
be concerned with dynamics on a phase space, TQ, and the dynamics will be 
those of a forced Lagrangian system. We will find that without active forcing, the 
system is dissipative and symmetric with respect to translations. This will allows 
us to perform Lagrangian reduction [3] to obtain an asymptotically stable point 
in TQ/M.. By perturbing this stable point we will be able to assert the existence 
of a limit cycle in TQ/K. using the "persistence theorem" of jl] and [9]. This 
construction can be seen as a terrestrial counterpart to recent work done by the 
author on swimming wherein the phase space for fluid structure interaction , A, 
exhibits an SE(3) symmetry which allows us to interpret swimming as a limit cycle 
on A/SE{3) [ID]. 

1.3. Acknowledgements. The notion of realizing the no-slip condition as a limit 
of viscous friction was brought to my attention by Dmitry V. Zenkov in a con- 
versation about the long time behavior of such systems. I was educated by Sam 
Burden on the role of limit cycles in model reduction for hybrid systems. A con- 
versation which involved banging our fists on tables over lunch convinced me that 
the role of phase space contraction was important for the passive control of these 
systems. Simultaneously Erica J. Kim asked how to view general locomotion as a 
limit cycle, and spurred my own pursuit on a project for viewing swimming as a 
limit cycle. The present paper can be seen as a terrestrial counterpart to this. I 
would like to thank Jaap Eldering for pointing out flaws in an earlier formulation 
of the model and providing much encouragement and interest in the project. Fi- 
nally, the initial stimulus to write this paper was given by Jair Koiller, who has 
been very supportive of my recent forays into biomechanics. This research was 
supported by the European Research Council Advanced Grant 267382 FCCA and 
NSF grant CCF-1011944. 

2. The model 

The model can be broken into two distinct components: the walker and the 
environment. We will discuss the model of the walker in a vacuum before we 
elaborate on how to model the ground, gravity, and the no-slip condition. 

2.1. A model of a walker (in a vacuum). The walker consists of three point 
particles of unit mass all connected by springs of unit stiffness with light viscous 
damping. These point particles move through IR 2 with positions q^ = (xi,yi) 
and velocities q = (xi,yi) for i = 1,2,3. If we let Q = IR 6 with coordinates 
q = (qi, q 2 , (fe) = yi, x 2 , y 2 , X3, 2/3), we can describe the walker in a vacuum as a 
Lagrange-mechanical system. Here kinetic energy is K(q, q) = 7,\\q\\ 2 and potential 
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energy from the springs is 

Uij(qi,qj) = r (yj(xi- Xj) 2 + (y t - y 3 ) 2 - i^t) 



where Uij(qi,qj) is the potential energy stored in the spring which connects mass 
% with mass j for some time dependent length £ij(t). The functions iij = 
do not control the shape of the walker directly, but instead determine an energy- 
minimizing shape and therefore control it indirectly 





Figure 2. Depicted is a cartoon of our walker. The grey cloud in 
the background represents the potential energy due to gravity plus 



the potential energy from the ground, U np of equation (2.2). 



The Lagrangian is given by 

L(v,q;t) = h\q\\ 2 - ^Uijiq^qj) 



i<j 



and we define the viscous force on the spring between particle i and j by 

Fij = (ij - Xi)dxi + (yj - yi)dyi. 

The equations of motion are q = (Yl^j =1 F^J — fX/i<j dUijj , and explicitly given 
in terms of g; and qj as 

Qi = Yl (ihj - Qi\\ ~ hjil))";, + qj - ft) 

where qi = (xi, y^) and Ay is the unit vector which points from mass i to mass j. 
As this is ill-defined when q^ = qj, the equations are ill-defined on a set of measure 
zero. 



2.2. A regularized model of the ground. Now we will present a model for 
the ground on which the walker will walk. The surface of the ground occupies 
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the line {y = 0} in IR 2 . Ideally, the ground is unpenetrable and imposes a no-slip 
condition. This is mathematically represented by the constraints 



(1) 
(2) 



^ > 

Xi = if yi 







for i = 1, 2, 3 where equation Q is the no-penetration condition and equation ^ 
is the no- slip condition. Both the no-penetration and no-slip conditions present 
challenges of a singular nature as they are conditions which abruptly "turn on" at 
y = and are inactive otherwise. It is precisely this "on/off" character which we 
will regularize. 

Firstly, the no-penetration condition can be realized by taking a potential energy 
that is infinite for y < and finite for y > 0. We will model this non-smooth 
potential by approximating it with a smooth one. In particular, the hyperbolic 
tangent 

tanh(a;) 



e~ 2x + 1 

serves as a smooth approximation of a step function and allows us to define the 



potential energy [/, 



np 



by 



U np (qi) = -/inptanh ^ 
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which has the gradient 



dU np (qi 
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1 



dyi 



and prevents the ith particle from falling through the floor. The parameter cr np 
controls the width of the region in which this force is active while /z np amplifies 
this force. As cr np approaches 0, this force becomes concentrated around the line 
{y = 0} and becomes infinitely strong; thus our model approaches an exact model 
when cr np approaches zero and /i np remains constant. 

The no-slip condition is similar to the no-penetration condition in that it is 
only active at {y = 0}. However, unlike the no-penetration condition, the no-slip 
condition is not derivable from a potential energy but instead can be viewed as a 
limit of viscous friction In particular, consider the viscous force on the ith 
particle given by 



^ns(0 , i; q% 



- tanh ( — 

2 V (Jnc 



This force "turns on" in the region near {y < 0} and dampens motion of the 
particles in that area. In other words, F ns resists slipping but does not outlaw it. 
As before, the parameter er ns controls the width of the boundary where this force 
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is active while /i ns controls the amplitude of this force. When <r ns goes to and 
p ns goes to oo we arrive at a no slip condition. 

Similarly, we desire to dampen bouncing. This is done with the viscous force 



Fdb(Qi,Qi) = -HdbXi 



- tanh 


'Vi_\ 


I 


2 ' 




~ 2_ 



dyi. 



Finally, we incorporate gravity via the potential energy on each particle given by 
U gT (x,y) = y. This potential energy imposes the force —dU gr (x,y) = (0,-1) on 
each particle. 

2.3. The full model. Now that we have established the Lagrangian of the walker, 
as well as the environmental forces imposed on it, we may finally provide the 
equations of motion. These equations of motion are obtained by adding the viscous 
forces, F ns and Fdb, and the potential forces, dU np and dU gr , to the equations for 
the walker in a vacuum. Therefore, the Lagrange-d'Alembert equations are given 
by 

- dU np (qi) + F ns (qi, q t ) + F dh (q i: q t ) - tZC/ gr (g i ). 
We may express these equations in a more familiar form by defining the total force 
F (q, Q) = ^2 F ij(<ii, Qj) + ft) + F dh (qi, qi) 

ij i 

and the total potential energy 

U(q) = ^2 Uijfai, qj) + Y U np(Qi) + U gI (qi) 

i<j i 

to write the equations as 

q = F(q,q)-dU(q) 

3. Qualitative Theory 

3.1. Reduction by Translation Symmetry. Consider the left action of R on 
IR 2 given by translating the ^-coordinate of each particle. We will denote this 
action by p : IR x IR 2 — > IR 2 . In this section we will show that our system is 
invariant with respect to this Lie group action. In particular, the action on Q 
which we also denote by p : IR x Q — y Q is given by 

P(s, (qi, q 2 , <?s)) := (p(s, qi),p(s, q 2 ),p(s, q 3 )) 

This action on Q can be lifted to TQ by the tangent lift. In coordinates, this is 
given by the map 

p(s) ■ (q,q) i—y (p(s) ■ q,q). 



LIMIT CYCLE WALKING ON A REGULARIZED GROUND 



7 



It can be observed that L is invariant with respect to this action on TQ; this is 
expressed formally by the observation 

L (p( s ) ■ (?,?)) = L(q,q) 

for any (q, q) G TQ and sei Additionally, the force F is really a map from TQ 
to T*Q and is also R invariant in the sense that 

(F(p(s) • (q,v 1 )),p(s) • (g,u 2 )) = (F(q, Vl ), (q,v 2 )) 

for any (g, vi), (q, V2) G TQ and any s G R. Given that the equations of motion 
are written in terms of L and F, there must exist an equation of motion on the 
quotient space TQ/M. The main goal of this section is to find these equations of 
motion on TQ/M. 

Let us begin by noting that R acts by shifting the x coordinate and so Q/M = R 5 . 
We can coordinatize Q/M by dropping the x coordinates of the third mass (this is 
one choice of many coordinate systems). This gives us coordinates {x\, y\, x 2 , y~2, Vz) 
on Q/M. and induces the coordinates for TQ/M = R 11 given by (xi,yi, x 2 ,y2, $3, V\, Wi,v 2 , w 2 ). 
We may construct the principal bundle projection n : TQ — > TQ/M., which makes 
the assignment 

xi = xi- x 3 , fa = y u vi=vi- v 3 , Wi = Wi 

x 2 = x 2 - x 3 , yi = yi, v 2 = v 2 - v 3 , w 2 = w 2 

yi = yi, v 3 = v 3 , w 3 = w 3 

The equations of motion on TQ/M can be written explicitly in terms of these 
coordinates. We find that the R- invariant potential energies Uij, U np : Q — > M may 
be written as functions Uij, Ufj : Q/M — > R, or more explicitly as 

Un(q) = I [VW^W+W^W ~ h2(t)) 
U np {qi) = -/inptanh ( — ) . 

We may define the total potential energy U = (^2^ + ^=1 U np (q,j), which 

is simply the manifestation of U on the reduced space Q/M. 

Moreover, the forces are also R- invariant. That is, for each force / : TQ — > 
T*Q introduced in the previous section, there is a map / : TQ/M — > (TQ/M)* 
characterized by the property 

(/(*"(?, v ))> v)} = (f(q, v), (q, v)) V(g, v), (q, v) G TQ. 
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Moreover, we may write the forces in terms of our coordinates on TQ/M.. In 
particular, they are: 



^31 = -vi{dxi + dx 3 ) + (w 3 - 

# 13 = v l( ix 3 + (wx - w 3 )dy 3 

F32 = —v 2 {dx 2 + dx 3 ) + (w 3 - w 2 )dy 2 

p 23 = v 2 dx 3 + (w 2 - w 3 )dy 3 

F\ 2 = (vi - v 2 )(dx 2 + dx 3 ) + (wi - w 2 )dy 2 

F 2 \ = (v 2 - Vi)(dxi + dx 3 ) + (w 2 - wx)dyi 

Fw,\ = -Vns{vi + v 3 ) ( ^tanh ( — 



Ftm,2 — 



-/ins (^2 + v 3 ) (- tanh ( — 



F ns ,3 = -Unsh Qtanh (J^~J ~ \ 
F d b,i = -Hb,$i Qtanh (^~J ~ d Vi 




where dxi, dx 2 , dx 3 , dy%, dy 2 , dy 3 are dual vectors to the space spanned by the co- 
ordinates Vi,v 3 ,v 3 ,Wi,w 2 ,w 3 - The dual vector dx 3 is not the differential of the 
x 3 coordinate because there is no x 3 coordinate on the space Q/M. The sequence 
of symbols Ll dx 3 " represents a dual vector to the space spanned by the coordinate 
v 3 . That dx 3 is not the differential of a coordinate is related to the fact that v 3 
is not associated to a base point and therefore is not a tangent vector in a tra- 
ditional sense. This is because TQ/M. is not actually a tangent bundle, but is a 
generalization of a tangent bundle known as a Lie algebroid [13J. 
We find that the base coordinates on Q/M evolve by 

di\ A dyi 

dt 1 ' dt 1 
dx 2 _ A dy 2 _ A 

dt ~ dt 
dyz 

If we define the total force to be F, then the evolution of v = (vi,wi, v 2 , w 2 ,v 3 , w 3 ) 
is given by 

d 4 = F-dU. 
dt 
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3.2. Hyperbolic Stable Points. In this section we will use the energy function 
on TQ /R as a Lyapunov function to assert the existence of a hyperbolically stable 
point for the passive dynamics. In particular, consider the energy function 

E(q^) = \\\n 2 + u(q). 

We can use the equations of motion to calculate the time derivative of E along a 
trajectory 

dE dE dq dE dv ,a-/-n * * , ,Av-\ A/ * A/- -\ - 

= ~d^dt + ^Jt = dU (l)- v + v -(- dU (<l) + F ^ v )) = F(q,v)-v 

It is simple to observe that ^ is negative definite in v since for any (g, v) such 
that (?, u) = ir(q,v), we find that 



■ v = F(q,v) ■ (q,v) = -|MI* + 7",^ + ttf»Wj 

where 



IHI2 = IHI a + X)/»(wK + /db(wK 
i=i 



/ns(y) = /"nstanh 
/db(l/) = /idb tanh 



(7,. 



dE 

— < - 
dt ~ 



- ■ 

V°"db/ 

Noting that X/i<j + w i w j < \\v\\ 2 we see that 

Y f™(yi) v i + fdh{yi)wj 

which is negative definite in v and w as long as yi > for % = 1,2,3. One can verify 
that U has isolated minima in Q/M. when /j, np is sufficiently large. Moreover, these 
minima occur in the region > for i — 1, 2, 3. If we denote one of these minima 
by ?mim then g min is an asymptotically stable fixed point because the energy is 
decreasing monotonically until v — w — 0, at which point the state of the system 
must be at a critical point of If. 

3.3. Time periodic perturbations. Given a dynamical system x = f(x) on a 
manifold M with a hyperbolically stable fixed point Xq in M, one can embed the 
system into a time-periodic augmented phase space S 1 x M by using the vector field 
= (1, f(x)). Then the orbit 70 (t) = {(0, £0) : $ e is a hyperbolically 
stable limit cycle for the system on S 1 x M. Moreover, given a time-periodic 
perturbation / 1— > f + egg, the limit cycle will continue to exist but be transformed 
into another limit cycle 7 e [I], [9] (see also "The Averaging Theorem" in j8j) . A 
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cartoon of this idea is depicted in figure [3] We will abuse notation by referring to 
both the limit cycle on S 1 x M and the projected cycle on M by the symbol 7. 




A B c 

Figure 3. (A) A one dimensional system with a hyperbolic stable 
point. The embedding of the system in (A) into a time-periodic 
augmented phase space results in a limit cycle as shown in (B). 
Finally, a limit cycle will persist under a time-periodic perturbation 
as a diffeomorphism of the original limit cycle. This is drawn in (C). 

In the previous subsection, we found a hyperbolically stable fixed point on 
TQ/M. In this section, we will perturb this system by allowing the lengths £y 
to be time periodic with period T. If these oscillations are small, we can expect to 
observe a stable limit cycle, j(t), also with period T in the augmented phase space 
TQ/1SL. However, this limit cycle is only a cycle on TQ/M. In general, if T(t) is a 
trajectory in TQ such that 7(2) = 7r(T(t)), then it is generally not the case that 
r(t) is periodic. In particular, any periodic trajectory 7 C TQ/M of period T is 
reduced from a trajectory (q, q)(t) G TQ such that 

(q,v)(t + T) = p(Ax)-(q,v)(t) 
where Ax e M is obtained from the v component of j(t) via 

f T 

Ax = / V3(r)dr. 
Jo 

The above integral may be viewed as a phase reconstruction formula with respect 
to the reduction by K symmetry. This relationship between trajectories in TQ 
and loops in TQ/M. is related to walking in the sense that the periodicity of j(t) G 
TQ/M. implies that (q,v)(t) G TQ satisfies 

(q, v)(t + kT) = p(kAx) ■ (q, v)(t). 

Simply put, the motion in TQ is relatively periodic in that each period is related 
to the previous period by a shift to the right by an amount Ax. For this reason, a 
trajectory (q, v)(t) G TQ is called a relative limit cycle if when j(t) = 7r((g, v)(t)) G 
TQ/M is a limit cycle. Moreover, this behavior in TQ is stable because the limit 
cycle j(t) G TQ/M is stable. Such behavior could resemble walking which is 
characterized by being relatively periodic and stable. 
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3.4. Early Numerical Results. I hypothesize that the shift Ax is generally non- 
zero. We can provide evidence in support of this hypothesis through the help of 
numerics. In this section, we use the stiff ODE solver of the python scientific 
computing library with an absolute tolerance of le-10 to integrate the equations 
of motion. With the parameters 

/i ns ^ 10.0 <T ns 0.1 

/idb | — > 1-00 <7db ^ o.i 

yU np (-)■ 10.0 (7 n p I — y 0.1 

along with the time-dependent lengths 

Z 13 = 1-0.3 sin (u (t-lj) 

Zi 3 = 1 — 0.3 sin (ut) 
Z23 = 3 — lu — Z13? 

where the frequency u = 2tt. This particular choice of parameters yields behavior 
which suggests that we are observing a relative limit cycle. A typical trajectory is 
depicted in figure [4] and the motion of the masses in a reference frame attached to 
the center of the body is depicted in figure |5j 

We observe an initial transient period followed by orderly periodic behavior. 
Over large times, this results in a constant drift of the walker towards decreasing 
values of x, while the y-coordinate of the walker stays just above y = 0. Figure 
[4] also provides a zoomed in view to illustrate what the evolution looks like on 
smaller time scales. On these smaller time scales we can see that the x-coordinate 
oscillates with an angular frequency of 2n and and a gradual but steady drift 
towards decreasing x, while the y-coordinate tends to oscillate with an angular 
frequency of 2tt without any drift at all. To further verify this behavior, we have 
also provided a plot of the trajectories in a moving reference frame such that the 
average x coordinate of the masses is in figure [5j One can interpret figure [5] 
as a plot of the trajectories of the masses in Q/R. One can observe some initial 
transient dynamics before these masses settle into cyclic behavior. The long-time 
cyclic behavior manifests as three thick circles in figure [5] located roughly at the 
coordinates (-0.5,0.3), (0.0,1.0), and (0.5,0.3). 

Altering parameter values and the initial conditions always produces behavior 
indicative of a relative limit cycle. However, the numerical results should be taken 
with a grain of salt, because multiple hyperbolic stable equilibria exist in Q/M.. 
In particular, given an equilibrium in Q/M one can permute the masses with one 
another to achieve another equilibrium. For example, the evolution depicted in 
figure [4] suggests that the system almost settles near one of these equilibrium 
configurations before the yellow mass and the red mass appear to swap positions 
between t = 10 and t = 20. One could imagine varying the initial condition only 
slightly, in such a way that yellow and the red mass do not sway and the system 
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20 40 60 80 100 

time 



Figure 4. This plot depicts the evolution of a typical trajectory. 
The solid lines correspond to the x coordinates of the three masses, 
while the dotted lines correspond to the y coordinates. We can 
clearly see a constant drift of the x-coordinate of each mass. If we 
blow-up a region we can directly observe behavior which suggests 
which suggests the system is captured in a stable periodic orbit in 
TQ/R. 



settles to a different relative limit cycle. This is just one particular example. 
Each equilibrium can be periodically perturbed to produce a limit cycle, and it 
is generally not the case that the same walking behavior will result for different 
limit cycles. A natural extension of this paper is understanding the consequences 
of having multiple hyperbolically stable equilibria. 

Secondly, while the trajectories appear to have converged to a limit cycle, it 
is not clear that absolutely all of the transient dynamics have left the system. 
There may be some slowly varying dynamics which have not dissipated. Further 
analytics may confirm or reject this posibility. Nonetheless, it is apparent that 
limit cycle-like behavior dominates fairly quickly in these systems. 



4. Conclusion & Future Work 

In this paper we have presented a model of walking where the walker consists of 
three point particles in IR 2 connected by springs, and the ground consists of the sum 
of a smooth potential energy and a smooth viscous force on the phase space, TQ. 
We then found that the system exhibits a symmetry with respect to translation of 
the x-coordinate. The reduced equations of motion on TQ/M. were dissipative, and 
lead to hyperbolically stable equilibria. By periodically perturbing the system we 
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Figure 5. This plot depicts the location of the three masses in a 
reference frame where the average x-coordinate of the masses has 
been translated to the origin. Therefore this plot depicts the trajec- 
tory of the masses in Q/M. Here we can observe limit cycle behavior 
in that each of the masses settles to stable oscillatory behavior after 
an initial period of transient dynamics. 

obtained limit cycles in TQ/M. in the vicinity of the equilibria. When these limit 
cycles in TQ/M. were lifted to relative limit-cycles in TQ we observed dynamics 
which resembled walking. Numerical simulations then provided further evidence 
for this interpretation. The work done so far suggests a number of avenues to 
pursue. In particular: 

(1) While we observe a limit cycle for the regularized model, it is not clear 
that this limit cycle persists in the limit as cr ns , <r np — > 0. Proving that 
this is (or is not) the case can be approached by viewing this as a singular 
perturbation problem. 

(2) As previously mentioned, while the numerics may make one initially happy, 
it is not clear that all of the transient dynamics have left the system. The 
influence of time-scales on this system is particularly important as the 
problem is likely exacerbated as the width of the regularization shrinks to 
zero. Deriving lower bounds on the decay rates for the transient dynamics 
could help resolve this issue. 

(3) While the limit cycle in the paper is stable, the size of the stability basin is 
not mentioned. Having a large stability basin is one method of achieving 
robustness and so having a lower bound for the radius if this basin would 
be useful. 
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In summary, it appears that regularized models are capable of exhibiting be- 
havior which resembles walking. Such models are open to classical techniques in 
dynamical systems, and allow one to view walking as a limit cycle in a reduced 
space in which case the absolute motion manifests as a phase shift. We hope that 
the simplicity of this characterization encourages further research down this av- 
enue which can enhance observations as well as complement the hybrid systems 
approach. 
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